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Selected results of a classical simulation of A*' bodies in strong interaction are presented. The static 
properties of such classical systems are qualitatively similar to the known properties of atomic nuclei. 
The simulations of collisions show that all observed reaction mechanisms in nucleus-nucleus collisions 
are present in this numerical simulation. The first studies of such collisions are in qualitative 
agreement with experimental observations. This simulation could shed a new light on energy deposit 
in heavy ion collisions: the excitation energy of each cluster is found lower than the energy of the 
less bound body in the cluster. Finally, the similarities of fragmentation pattern with those obtained 
with a statistical code indicate that an unambiguous link can be established between the statistical 
and the dynamical descriptions of reaction mechanisms. 
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During the last decades, heavy ions collisions have been 
widely used and intensively studied experimentally to de- 
termine the equation of state of nuclear matter. Many 
reaction mechanisms have been identified and a lot of 
theoretical works has been induced. Schematically, the 
incident energy domain has been divided in three parts: 
the low energy range, up to 10'--^20 A.MeV, where fusion 
and deep inelastic processes are dominant; the high en- 
ergy range, above 100 A.MeV, where the nucleonic and 
sub-nucleonic degrees of freedom start to play a signif- 
icant role; and the intermediate energy range, from 20 
A.MeV to 100 A.MeV, where binary processes and multi- 
fragmentation take place. At each energy range, the reac- 
tion is schematically described as a two steps process: an 
entrance channel during which excited nuclei are formed, 
and a decay stage where the excited fragments cool down 
by particle or fragment emissions. 

Theoretically, the entrance channel of the reaction is 
mainly described by dynamical models. At low energy, 
the attractive character of the interaction is dominant 
and the entrance channel is described mainly by mean- 
field approaches ||l|. At high energies, the repulsive part 
of the interaction and its nucleon-nucleon character are 
dominant and cascade models are mainly used to de- 
scribe the first moments of the collision At inter- 
mediate energies, the attractive part and the repulsive 
part of the interaction interfere. Transport models are 
widely used to describe the entrance channel at this en- 
ergy range ||, |^, ||, ^, 0, ^, ^ • They contain a mean-field 
part through the one-body evolution of the system, and 
a nucleon-nucleon part through the collision term. At 
all energy ranges, the decay stage is mainly described 
by statistical decay models (GEMINI SIMON 0, 
SMM (l|, MMMC |l|, QSM @, EES]l|, etc.) which 
assume that the nuclei formed during the first dynamical 
step are equilibrated. 

With this scheme, it is very hard to have a global 
and consistent description of nucleus-nucleus collisions. 
Although the nucleon-nucleon interaction is the same. 



whatever the incident energy and whatever the reaction 
time, different models are used depending on the energy 
range and the reaction time. The links between these 
different models are not obvious. For example how the 
liquid-drop parametrisation used in the statistical decay 
code can be deduced from the parameters of the interac- 
tion used in a transport code? What is the link between 
the incompressibility modulus of nuclear matter K^o and 
the nucleon-nucleon cross section cr„„? Is the statistical 
description still justified when the thermalisation time is 
close to or longer than the life time of excited nuclei? 

The present status results mainly from the impossi- 
bility to solve the full nuclear many-body problem. To 
solve it, approximations are made, depending on the en- 
ergy range and on the reaction time. Even after simpli- 
fications, the resulting equations do not have analytical 
solutions and they are often solved numerically. It seems 
hard nowadays to connect these different models to have 
a global description of the processes involved in nucleus- 
nucleus collisions. 

Nevertheless, some attempts were made to describe 
globally nucleus-nucleus reactions at intermediate ener- 
gies by using (semi-) classical molecular dynamics codes 
0, 1^, |18|, |l|, |0[ |l|, H. The main advantage of 
these codes is that the classical many-body problem 
can be solved numerically with a high accuracy with- 
out any assumption. The main drawback is that the 
quantum character of the system is ignored. A lot of 
work has already been done with the help of such codes 
to study the mechanism of fragment formation in multi- 
fragmentation. Could these codes do better by describ- 
ing the whole reaction process, from the very beginning 
of the reaction up to the decay stage? Are they able to 
make a link between the different energy ranges? Which 
processes are accessible with the simplest hypothesis? 

This paper will show that classical molecular dynam- 
ics simulations are interesting tools to establish the links 
between the different approaches used in nucleus-nucleus 
collision. In the first section, the Classical N-Body Dy- 
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namics code will be briefly described. The static prop- 
erties of stable N-body systenis will be studied in the 
second section. Some examples of reaction mechanisms 
and some analyses of cluster-cluster collisions are given 
in the third section. Finally, conclusions will be drawn. 
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FIG. 1: Shape of the interaction used in the CNBD simulation 
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on a Compaq DS20 computer under the UNIX True64 
operating system ranges from « 30 to « 400 seconds, 
depending mainly on the impact parameter. The energy 
gap between the beginning and the the ending simulation 
time is lower than 0.001%. This simulation has five free 
parameters: four linked to the physics (the interaction) 
and one linked to the numerical algorithm (e). 

Since one wants to study the simplest case, neither long 
range repulsive interaction nor quantum corrections like 
a Pauli potential have been introduced [Q. Addition- 
ally, no statistical decay code is applied on the excited 
fragments formed during the collision. The final prod- 
ucts have to be regarded as "primary" products which 
will decay afterwards. 

In order to avoid any confusion with nuclear physics, 
the units used here are called Simulation Units and 
noted S.U.. The distance will then be in Distance Sim- 
ulation Units {D.S.U.), the energies in Energy Simula- 
tion Units {E.S.U.), the velocities in Velocity Simulation 
Units (V.S.U.) and the reaction time in Time Simulation 
Unit {T.S.U.). We will be only interested in the rela- 
tive evolutions of the observables and in their link to the 
properties of the stable systems. The main goal of the 
present work is not to reproduce the experimental data 
of nucleus- nucleus collisions, but rather to see to which 
extent this simple simulation is qualitatively similar, or 
not, to experimental data. 



I. DESCRIPTION OF THE CODE. 



II. STATIC PROPERTIES OF "GROUND 
STATES". 



Let us start by describing the Classical N-Body Dy- 
namics code (labeled CNBD) used in this article. The 
basic ingredients of such a code are very simple. The dy- 
namical evolution of each body of the system is driven by 
the classical Newtonian equations of motion. The two- 
body potential used in the present work is a third degree 
polynomial whose derivatives are null at the range ri and 
at the distance of maximum depth tq. The depth value 
is Vmin and the value at r = is finite and equal to 
Vq. The shape of the two-body potential is shown on 
figure |[ This potential has the basic properties of the 
Lennard- Jones potential used in other works 22 : a 
finite range attractive part and a repulsive short range 
part. To follow the dynamical evolution of the system 
an adaptative stepsize fourth-order Runge-Kutta algo- 
rithm is used [ p3[ . The main difference with other works 
is that the time step At can vary: if the potential varies 
strongly. At is small and when the potential varies gently. 
At becomes larger. This allows a very high accuracy with 
shorter CPU time than for fixed time step algorithms. It 
requires an additional simulation parameter e which is 
adjusted to ensure the verification of conservation laws 
(energy, momentum, angular momentum) with a reason- 
able simulation time. For an e value of 10"^, typical CPU 
times for a collision of two clusters with 50 particles each 
with an ending time equal to 200 time simulations units 



Once the basic ingredients are defined, one can build 
stable systems. Since the two-body potential only de- 
pends on the distance between the two bodies, such sys- 
tems are small crystals. The "ground states" of such sys- 
tems are defined as the configuration in position space 
which minimizes their total energy. This is obtained by 
using a Metropolis simulated annealing method |2^ . The 
locations of particles obtained this way are very close to 
those obtained by using a basin-hopping algorithm for 
Lennard- Jones clusters ^ . 

On top-left of figure | is displayed the energy per body 
EBind/N of these "ground states" as a function of the 
number N of bodies in the cluster (solid line). The 
dashed and dotted line corresponds to a fit using a liquid- 
drop formula. The dashed line corresponds to the energy 
of the less bound body ELessBound m the cluster. It will 
be seen that this energy seems to play a particular role 
in cluster-cluster collisions. The overall dependence of 
Esind/N with N is very similar to what is seen in nu- 
clear physics. The main difference is seen for high values 
of N because of the absence of a Coulomb-like interac- 
tion: Esind/N continues to decrease with N whereas it 
increases for nuclei. These values are close to those found 
by using the basin- hopping algorithm . 

The dependence with iV of the root mean radius radius 
r{N) = V < r"^ > of N-body clusters is shown on the 
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FIG. 2: Energy per body (upper left panel) and root mean square radius (upper right panel) of N-body clusters as functions of 
N . For these two panels, the full line correspond to the obtained values and the dashed and dotted line to a fit. On the left 
panel, the dashed line corresponds to the energy of the less bound body in the cluster. Lower row: zero temperature "equation 
of state" for different N values (left) and "incompressibility modulus" as a function of N (right). 



top-right panel of figure || (solid line). The dashed and 
dotted line corresponds to a fit using a TVs term. Here 
again, this is similar to what is known for nuclei. The 
main difference is the rg term which is due to the large 
size of the "repulsive core" compared to the range of the 
attractive part of the interaction. 

Since a size can be defined, an effective density can 
be calculated. One can stretch and squeeze the N-body 
cluster and build an effective "zero temperature equa- 
tion of state" . This is displayed on the low-left panel 
on figure || for various system sizes. The density pe// 
on this abscissa is the density relative to the density for 
the "ground state". The curvature of this "equation of 
state" can be computed and defined as the "compressibil- 
ity modulus" K{N) of the system. Its variations with N 
are displayed on the low-right panel of figure 0. One can 



see on these two plots that the "equation of state" and 
K{N) are strongly dependent on N for values of N below 
20^30 and then less dependent on N above. One could 
be able to define an "equation of state of infinite mat- 
ter" by computing the limit of these evolutions for huge 
values of N. Here again, the behavior of these classical 
N-body systems is very similar to that of nuclei. 

This first study shows that the classical N-body clus- 
ters have strong similarities with atomic nuclei. This 
suggests that all the parametrisation used to describe nu- 
clei (liquid drop parametrisation, radii, equation of state, 
etc..) could be directly deduced from the parameters of 
the two-body interaction. Since this interaction remains 
unchanged, the differences seen for different clusters has 
to be attributed only to N and to the geometrical con- 
figuration of the bodies in the clusters. 
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FIG. 3: Pictures of cluster- cluster collisions for different reaction mechanisms. On each panel, Nproj is the projectile size, 
Ntarg the target size, V the relative velocity in V.S.U. between the target and the projectile , b the impact parameter in D.S.U. 
and t the collision time in T.S.U.. 



III. CLUSTER-CLUSTER COLLISIONS. 



After studying the static properties of stable N-body 
systems, collisions between clusters can be investigated 
to see what kind of reaction can be obtained. Roughly 
20,000 collisions have been generated for different system 
sizes, different entrance channel asymmetries and differ- 
ent impact parameters. For each collision, the orientation 
of the inertia axis of clusters are randomly chosen. This 
avoid to run twice the same collision if the impact param- 
eter, the target and the projectile size and the incident 
energy are the same for two different collisions. The sim- 
ulations have been performed for the following systems: 
13 + 13, 34 + 34, 50 + 50, 100 + fOO and 18 + 50. 
The incident energies have been chosen in a way that the 
available energy in the center of mass frame ranges from 
an energy far below the binding energy of the fused sys- 
tem (« 30 E.S.U) to an energy well above the binding 



energy of the fused system (sa 120 E.S.U.). As it will be 
shown, these collisions can be studied in the same way as 
the experimental data of nucleus-nucleus collisions are. 
These collisions can be seen as "numerical experiments" . 
Only few examples of such analyses will be shown. More 
detailed analyses will be done in forthcoming articles. 
Firstly, a small list of reaction mechanisms obtained will 
be done. Secondly, the excitation energy stored in excited 
clusters will be studied and a link to the static proper- 
ties will be done. Finally, the question of the statistical 
description of this dynamical model will be briefly ad- 
dressed. 



A. Reaction mechanisms. 

On figure || are displayed the pictures of cluster-cluster 
collisions obtained at a fixed time. The ending time 
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FIG. 4: Wilczinky plots for N = ?,A + N = 34 collisions at different available energies in the center of mass (see text). 



of these collisions is t — 200 T.S.U.. On each picture, 
the bodies which were originally belonging to the projec- 
tile are in dark grey and those which were belonging to 
the target are in light grey. The most striking observa- 
tion is that all kinds of reaction mechanisms observed in 
nucleus-nucleus collisions seem to be present in this very 
simple simulation. One can see low energy processes like 
the fusion/evaporation process (first row left panel), pure 
binary collision (second row left panel), stripping/pick- 
up mechanism (second row right panel) and deep inelas- 
tic collisions (third row right panel). Intermediate en- 
ergy processes, like neck formation and break-up (first 
row right panel) and multi-fragmentation (third row left 
panel), are seen. Finally, high energy processes like the 
participant/spectator process (lower-most panel) are also 
seen. 

This similarity can also be seen on the so-called 
Wilczynski plots shown on figure |[ These plots dis- 
play the correlation between the flow angle 6 flow and the 



total kinetic energy of the clusters Ekin = where 
El is the kinetic energy of cluster i in the center of mass 
frame. Two bodies are assumed to belong to the same 
cluster if they are in interaction, i.e. if their relative dis- 
tance is below the range ri of the interaction (see figure 
|l|). This algorithm of cluster recognition is the simplest 
one and is called Minimum Spanning Tree algorithm in 
other works ^ . The most dissipative collisions cor- 
respond to the smallest Ekin values: the available energy 
is converted in internal energy of clusters. These plots 
were built for a = 34 projectile colliding a = 34 
target at four available energies Ecm/N in the center 
of mass: E„„/N = 60 E.S.U., E„„/N = 90 E.S.U., 
Ecm/N = 120 E.S.U. and Ec^/N = 30 E.S.U. which 
corresponds respectively to the energy of the less bound 
body for the fused = 68 cluster, to the binding energy 
per body for = 68, to the energy of the most bound 
body for A^ = 68 and to an energy below the energy 
of the less bound body for A^ = 68. For each energy, 
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FIG. 5: Excitation energy of clusters versus their parallel velocity for N — 34 + N — 34 collisions at different available energies 
in the center of mass. On each panel, the full line corresponds to the expected evolution for a pure binary process, the dashed 
horizontal line corresponds to the energy of the less bound body for N = 68 and the circle corresponds to the expected values for 
the fused system. 



1,000 collisions have been computed assuming a flat im- 
pact parameter distribution ranging from b — D.S.U. 
to the sum of the two cluster radii plus the range of the 
interaction {b^ax = r{Nproj) + r{Ntarg) + n, where ri is 
the range of the two-body interaction). In the analyses, 
each collision is weighted assuming a triangular impact 
parameter distribution between and bmax {weight ex. b) . 

At Ecm/N = mE.S.U. (upper left panel), the flow 
angle is always negative. This means that the projectile 
and the target like clusters are deflected to the opposite 
direction relative to their original one. At this energy, 
the attractive part of the interaction is dominant. For 
the lowest Ekin values, all possible values of 9 flow are 
covered: this corresponds to the fusion/evaporation pro- 
cess. At Ecm/N = QQE.S.U. (upper right panel), the pic- 
ture is changed and the range of 9 flow values is smaller 



than for the previous energy: the repulsive part starts to 
act. The fusion/evaporation area (small E^in values, all 
9 flow values) is still present. For the two highest energies 
(lower panels), the picture is roughly the same. For the 
less dissipative collisions, the attractive part is still dom- 
inant (negative 9 flow values). But when the dissipation 
increases, the repulsive part becomes dominant and the 
projectile and the target bounce on each other (positive 
9 flow values). For the most dissipative collisions, only 
positive 9 flow values are seen which indicates the disap- 
pearance of the fusion/evaporation process. These evolu- 
tions are qualitatively similar to what is seen in nucleus- 
nucleus collisions (see for example 28 ). 



These two studies indicate that the description of the 
reaction mechanisms in terms of mean-field at low ener- 
gies and in terms of nucleon-nucleon collision at high en- 
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FIG. 6: Excitation energy of clusters versus N for N = 34 + N — 34 collisions at different available energies in the center of 
mass. On each panel, the full line corresponds to the energy of the less bound body, the dashed line corresponds to the binding 
energy per body and the circle corresponds to the expected values for the fused system. 



ergies can be deduced from the properties of the nucleon- 
nucleon interaction only. The similarities of this classical 
N-body simulation with nucleus-nucleus collisions sug- 
gest that, apart for the quantum mechanics effects, an 
effective description of nucleus- nucleus collisions could be 
obtained with a reduced number of parameters. Provid- 
ing the interaction has a finite range attractive part and 
a short range repulsive part, the overall behavior of the 
N-body systems seems to be independent of the values of 
the parameters of this interaction. For a more quantita- 
tive agreement with nucleus-nucleus collisions, additional 
physical ingredients (quantum mechanics. Coulomb in- 
teraction) have to be included. 



B. Energy deposit in clusters. 



Let us now focus our study of cluster-cluster collision 
on selected topics. Of particular interest is the excita- 
tion energy stored in clusters. How this energy can be 
linked to the available energy and how it is linked to the 
properties of the "ground state" characteristics of these 
clusters? 

One can plot for example the evolution of the excita- 
tion energy E* /N of the cluster with its parallel velocity 
V7/ for iV = 34 + TV 34 collisions (figure |) for the 
whole impact parameter range. The excitation energy 
of each cluster is simply the difference between the total 
energy (potential plus kinetic) and the "ground state" 
energy of the cluster. It reads: 
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N 

E*=J2 ^k^n + E ^('^^^■) - EB^nd{N) (1) 
2—1 i^jyi 

where E].^^ is the kinetic energy of the body i in the 
cluster center of mass frame, rij the relative distance be- 
tween the bodies i and j, V{rij) the potential energy 
and EBind{N) the binding energy of the "ground-state" 
of the cluster with N bodies. This excitation energy is 
determined at the end of the calculation corresponding 
to t = 200 T.S.U.. This energy is very close to the one 
obtained at the separation time of the clusters (the small- 
est time at which clusters can be identified) , since in this 
time range the evaporation is very weak and the clus- 
ters have no time to cool down significantly On 
each panel of figure S, the full line corresponds to the 



expected correlation between E*/N and V/ / for a pure 
binary scenario (the excitation energy is only due to the 
velocity damping of each partner) , the horizontal dashed 
line to the energy of the less bound body for the fused 
system = 68 and the small circle is centered around 
the expected values of velocity and excitation energy for 
the fused system. 

At Ecm/N = 30 E.S.U., the points are slightly be- 
low the full line. This means that the excitation en- 
ergy is strongly linked to the velocity damping. The 
small shift is due to mass transfers between the projec- 
tile and the target, and to promptly emitted clusters. 
The area corresponding to the fused system is well popu- 
lated showing that a complete fusion process occurs. At 
Ecm/N = 60E.S.U., the distribution of points is roughly 
compatible with the pure binary process hypothesis, but 
the complete fusion process area is empty. The horizontal 
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FIG. 8: Dalitz plots obtained with SMM calculations for different excitation energies. The source has a charge Z 
A = 205 and the freeze-out density is pf.o. = po/3. 



line seems to be an upper limit to the excitation energy 
which can be stored in clusters. Clusters with rather 
small excitation energies are found at velocities around 
the center of mass velocity. For the two highest energies, 
this trend is enhanced: a set of clusters are around the 
full line, and when the full line is above the horizontal 
line, one finds clusters at small excitation energies around 
the center of mass velocity. The energy of the less bound 
body seems to be a limit to the excitation energy which 
can be stored in these clusters. 

This can be more clearly seen when the excitation en- 
ergy E* /N is plotted as a function of A^, as in figure ||. 
On each panel, the full line corresponds to the energy of 
the less bound body ELessBound m the cluster and the 
dashed fine to the binding energy per body E^nd/N. As 
in figure |^, the small circle corresponds to the expected 
values for the fused system. At Ecm/N = 30 E.S.U., 



the area corresponding to complete fusion is filled and 
all the available energy can be stored as excitation en- 
ergy. But for higher energies, one can clearly see that for 
each fragment size, E* /N never overcomes ELessBound- 
At Ecm/N = 60 E.S.U. clusters with sizes higher than 
the projectile size and the target size can be seen. This 
area corresponds to an incomplete fusion process. For the 
two highest energies (Ecm/N — 90 E.S.U. and Ecm/N = 
120 E.S.U.) , the plots are almost identical: there is no 
more fusion and the clusters are smaller than the target 
and than the projectile. One can notice that E* /N never 
reaches the binding energy E^ind/N except for small clus- 
ters where Eund/N and ELessBound are equal. 

This limitation of excitation energy can be understood 
quite easily. The less bound body remains bound to the 
cluster only if its total energy is negative, i.e. its kinetic 
energy due to the excitation is below its potential one. If 
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one assumes that the excitation energy is roughly equally 
shared over all bodies in the cluster, when the kinetic 
energy balances the potential energy of the less bound 
body, this body is no more bound to the cluster and 
can escape. To be observed for a long time, the excited 
cluster must have an excitation energy per body below 
the energy of the less bound body. 

The mechanism of energy deposit in classical N-body 
clusters seems to be the following one: the excitation 
seems to be mainly driven by the velocity damping of 
the two partners and to a lesser extent by exchanges of 
bodies between them. Once the energy of the less bound 
body is reached, unbound bodies and/or clusters escape 
quickly and keep an excitation energy per body below 
the energy of the less bound one. As a consequence, the 
highest energy deposit per body can only be obtained at 
energies close to ELessBound- For higher available ener- 
gies, the system fragments quickly, leaving rather "cold" 
clusters around the center of mass velocity. This could 
be an explanation to the quite low excitation energies of 
clusters found in central collisions for the Xe + Sn system 
at 50 A.MeV ||2^. This subject will be more completely 
covered in a forthcoming article. 



fission and multi-fragmentation when the excitation en- 
ergy increases. 

The strong qualitative similarity between these two 
pictures suggests that it could be possible to build a sta- 
tistical decay code which would give the same results as 
the dynamical simulation. Of particular interest would 
be the study of the relations between the parameters of 
the statistical decay code (size of the source, excitation 
energy, deformation, radial flow, freeze-out volume, ex- 
ternal constraints) with the parameters of the dynamical 
one (projectile and target sizes, impact parameter, pa- 
rameters of the two-body interaction). This could also 
allow to check under which conditions the statistical de- 
cay code can be applied. This study could finally allow 
us to have a consistent description of cluster-cluster colli- 
sions, where the parameters of the statistical decay model 
are deduced from the parameters of the dynamical simu- 
lation. In this case, the cluster-cluster collisions could be 
described completely with a reduced number of parame- 
ters. 



IV. CONCLUSIONS 



C. statistical description of collisions. 

Let us end with fragment size analyses of cluster- 
cluster collisions. In nucleus-nucleus collisions studies, 
such analyses are very often used to fix the parame- 
ters of statistical models and to verify the compatibility 
of statistical decay models with experimental data (see 
for example [^0|, |l], ^). The so-called Dalitz plot for 
central (0 < b/hmax < 0-1) = 34 -I- A^ = 34 col- 
lisions are shown on figure Each panel corresponds 
to a fixed available energy in the center of mass rang- 
ing from Ecm/N = 30 to 120 E.S.U.. For this analysis, 
only the three heaviest clusters are taken into account. 
Each event is associated with a point in this plot. The 
distance of one point with respect to each edge of the 
triangle is proportional to the size of each of the three 
clusters. The corners correspond to events with a large 
size cluster and two small ones (fusion/evaporation pro- 
cess), an event in the middle of an edge corresponds 
to an event with two equal size clusters and a small 
one (fission process or binary collision) and the center 
of the triangle corresponds to three equal size clusters 
(multi- fragmentation process). One can see that when 
the available energy increases, the reaction mechanisms 
goes continuously from fusion/evaporation to binary col- 
lisions and finally to multi-fragmentation and/or vapor- 
ization. 

This picture is qualitatively very similar to the one ob- 
tained in SMM calculations Q , describing the decay of 
a single source, as shown on figure ^. In this case, each 
panel corresponds to an excitation energy. As in figure ^, 
the lowest energies correspond to fusion/evaporation like 
processes, and goes with a continuous transition towards 



As it has been seen in this brief overview of the charac- 
teristics and the reaction mechanisms of classical A^ bod- 
ies in strong interaction, there are strong qualitative simi- 
larities between these systems and the atomic nuclei. The 
static properties and the reaction mechanisms observed 
for the atomic nuclei and for these classical clusters are 
found very close to each other. This could mean that 
the experimental observations made for nucleus-nucleus 
collisions are mainly governed by the N-body character 
of the system and by the overall shape of the two-body 
interaction (finite range attractive part and short range 
repulsive part). More physical ingredients (Coulomb in- 
teraction, quantum mechanics) are of course necessary if 
one wants to have a quantitative agreement with exper- 
imental data. But it is surprising to have such a quali- 
tative agreement while essential physical ingredients are 
missing in the simulation. 

Since all kind of reaction mechanisms are observed 
in these classical simulations, from the low energy fu- 
sion/evaporation processes to the high energy partici- 
pant/spectator processes, they may also allow to con- 
nect in a consistent way the mean-field and the nucleon- 
nucleon approaches. One could for example study the 
relation between the incompressibility modulus and the 
size of the repulsive part of the interaction {K^o and (t„„ 
in the nuclear case). It is well known in transport cal- 
culations that one can "stiffen" the effective equation of 
state by just increasing the value of the nucleon-nuclcon 
cross section (7„„ Q. 

Finally, such simulations may reconcile two approaches 
which were up to now opposed in nuclear physics: the dy- 
namical description and the statistical description of frag- 
ment emissions. The dynamical models are unique tools 
to establish the link between the parameters of the inter- 
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action and those of the statistical models. Additionally, 
one can determine under which conditions the statistical 
models can be used, and what arc the true meanings of 
the parameters used in statistical models (temperature, 
chemical potentials, freeze-out volumes, etc.). This study 
may allow to have a complete and consistent description 
of colliding N-body systems with a reduced number of 
free parameters. This could be not only useful for nu- 
clear physics, but also for other fields of physics like for 
example the cluster physics. 
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